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A general formalism of X-ray scattering from different kinds of surface morphologies is described. 
Based on a description of the surface morphology at the atomic scale through the use of the paracrys- 
tal model and discrete distributions of distances, the scattered intensity by non-periodic surfaces 
is calculated over the whole reciprocal space. In one dimension, the scattered intensity by a vici- 
nal surface, the two-level model, the N-level model, the faceted surface and the rough surface are 
addressed. In two dimensions, the previous results are generalized to the kinked vicinal surface, 
the two-level vicinal surface and the step meandering on a vicinal surface. The concept of Crystal 
Truncation Rod is generalized considering also the truncation of a terrace by a step (yielding a 
Terrace Truncation Rod) and a step by a kink (yielding a Step Truncation Rod). 



INTRODUCTION 



The control of the morphology and structure of surfaces is one key point to elaborate high quality devices, since 
defects such as steps, kinks or more generally surface roughness may drastically influence their physical properties. 
Defects may also be used to promote catalytic reactions or as preferential nucleation sites for the self-assembling 
£--~. ■ of nanostructures [3, e. g. in view of ultra-high density recording. In this context, the development of surface 
characterization tools is highly appreciated. If scanning probe microscopy is the technique of choice to get local 
real space information, statistical averages and in situ experiments are advantageously performed with diffraction 
techniques. Numerous real time studies of thin films growth, coarsening or roughening have been performed, for 
instance with X-ray SSJHli Spot Profile Analysis Low Energy Electron Diffraction (SPA-LEED) Helium Atom 

Scattering (HAS) @, |{J, LLi| or Reflection High Energy Electron Diffraction (RHEED) [llj techniques. 
In particular X-ray scattering have the advantage over the other techniques of a weak scattering cross section which 
makes possible the quantitative analysis of the diffraction patterns in the framework of the kinematic approximation. 
The purpose of this paper is to revisit the scattering theory by disordered stepped and kinked surfaces in the context 
of the development of grazing incidence X-ray studies in surface science. It extends elder works which essentially dealt 
with Helium and electron diffraction. 

The main problem in the calculation of the diffraction by disordered surfaces is to have a statistical description 
of the morphology and to be able to calculate the scattered intensity. If the morphology is perfectly periodic, the 
scattered intensity is calculated in a straightforward way. The structure factor of the unit cell which repeats from 
the surface to the bulk, i.e. the Fourier transform of its electronic density, has to be evaluated. Then the scattered 
intensity is the modulus square of the sum of the amplitude scattered by all the unit cells. In the special case of 
surface diffraction, a semi-infinite sum has to be done in the direction perpendicular to the surface. If one considers 
the simple case of a cubic lattice, the intensity reads: 
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Where a, b and c are the lattice parameters along the coordinates x, y and z. q x , q y and q z are the reciprocal space 
coordinates and S is the Dirac distribution. 

This well known result tells us that the intensity is not zero in between Bragg peaks in the direction perpendicular 
to the surface. Such scattering rods are called Crystal Truncation Rods (CTR) [l2| and arise because the crystal 
is truncated. However, in general, the surface morphology deviates from this ideal perfectly ordered case. The 
introduction of disorder may lead to significant changes of the scattered intensity, particularly in anti-Bragg condition 
where the intensity is proportional of that scattered by half a lattice plane, and to an increase of calculation complexity. 
It is necessary to consider the scattered intensity by objects for which sizes and distances from one to another are 
correlated. For instance in the case of a vicinal surface, the distance separating two nearest neighbor terraces from 
their respective centers is completely determined by their widths. Therefore it is impossible to describe the scattered 
intensity as the multiplication of the square modulus of the form factor of the terrace width (averaged over the widths 
distribution) by an interference function. Several authors have addressed this problem and some analytical results 
have been obtained assuming generic distributions of terrace widths [3, EE EE EE EE EH- We propose here a 
method based on the paracrystal model ID, 2^, 21 1 and an appropriate choice of the elementary scattering object at 



the atomic scale to evaluate the scattered intensity from stepped or kinked surfaces. Previous results of HAS and 
electron diffraction, i.e. the two-level model, the N-level model, the vicinal surface and the randomly stepped surface 
are extended to X-ray scattering. In comparison with these elder results the contribution of deep layers gives rise 
to a CTR prefactor. This one-dimensional approach is further expanded to two dimensions for some specific surface 
morphologies. All calculations are performed assuming the atomic structure of matter (discrete values of the terrace 
size and the kink-kink distance) which allows evaluating the scattered intensity in the whole reciprocal space and not 
only around Bragg peaks. This point is of importance since surface structure characterization is usually performed 
measuring the CTR or during growth by measuring the scattered intensity in anti-Bragg conditions. Therefore Small 
and Wide-Angle scattering data may be analyzed with the same model. The concept of CTR is extended to the 
truncation of a terrace by a step (Terrace Truncation Rod) and to the truncation of a step edge by a kink (Step 
Truncation Rod) thanks to the terrace-step-kink model used for the description of the surface morphology. 

This paper is organized as follow. The method is first presented for a classical one-dimensional analysis of surfaces. 
Different surface morphologies are successively treated: (i) the vicinal surface, (ii) the two-level model, (iii) the N-level 
model and (iv) the randomly stepped surface (co-level model), illustrating the choice of the elementary scattering 
object as well as that of paracrystals. The cases of faceted surfaces are next treated, demonstrating the generality 
of the formalism. Follows a generalization to two-dimensions for some specific cases: (i) the kinked vicinal surface, 
(ii) the two- level vicinal surface and (iii) the co-level vicinal surface (or step meandering on a vicinal surface). In a 
last part, a method is briefly proposed to take into account in the calculation the intensity scattered by the elastic 
distortions induced by periodic steps and kinks below the surface. 

ONE DIMENSIONAL ANALYSIS OF SCATTERING FROM SURFACES 

The vicinal surface 



The usual approach for the calculation [22], |23| of the X-ray diffraction by a vicinal surface with constant step 
height is based on a description of the surface as a stacking of terraces along the x-axis (Fig [T]) . Then the scattered 
intensity is averaged over all possible configurations of the step positions. Despite its simplicity from a formal point of 
view, it is necessary to take into account the correlations between the terrace widths and the distances that separate 
them. However a precise analysis of vicinal surfaces shows that an elementary object without any size distribution 
can be extracted from the surface morphology. This object consists of a semi-infinite plane parallel to the terraces 
(along the x-axis), starting at the step edge (Fig ^ and penetrating into the bulk. Concerning the terrace width 
distribution and the average over all the configurations, the use of the auto-correlation function of the step positions 
in the framework of the paracrystal model is straightforward. Therefore this method allows separating in a simple 
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way the role played by the step structure from that played by the terrace width distribution. In the framework of 
this model, the scattered intensity is the product of the square modulus of the form factor of the elementary object 
F vicina i{q x ,q z ) by an interference function S(q x ,q z ). 

I(Qx,Qz) = \ F vicina i(q x ,q z ) | 2 xS(q x ,q z ) (2) 

The form factor of the elementary object is a semi-infinite sum, which yields the usual Crystal Truncation Rod 
behavior parallel to the terrace width: 
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where F ce u(q x ,q z ) is the structure factor of the unit cell, a is the lattice parameter along the x axis and /j, is the 
projected absorption coefficient along the x axis. One must remember that the above formula is not correct near Bragg 
peaks of perfect crystals due to multiple diffraction effects. The second part of the calculation concerns the interference 
function, i.e. the Fourier transform of the auto-correlation function of all the step positions. The autocorrelation 
function, g(x, z), is obtained as an infinite summation over n of the probabilities P n (x, z) of having a n th neighbor 
step (A n ) at a distance (x,z) from a step A arbitrary located at the origin (P (x,z) = S(x) x S(zj). 

oo 

g(x,z) = J2 p n( x ,z) (4) 

n— — oo 

To calculate each P n (x, z) term, we assume that the distance between nearest neighbor steps is described by a 
probability law, the terrace widths distribution Pr(x), and that there is no correlation between the sizes of neighboring 
terraces. Within this assumption, P n (x, z) can be obtained step by step building the vicinal surface. The probability 
for the first neighbor step A\ to be located at the (x,z) coordinate is P\{x, z) — Pt{x) x S(z — c), where c is the 
step height. The probability of having a second neighbor step A2 at (x',z') results from the multiplication of the 
probability of having a first neighbor step Ai located at (x,z) by the probability of having A 2 and A\ separated by 
the vector (x'-x,z'-z) integrated over all intermediate A\ positions. This probability is simply the convolution product 
of two probability laws : Pi (a;, z) ® Pi(x, z) — Pi(x, z)® 2 . This result can be generalized for the n th neighbor step : 

P n (x,z) = [Pi(x,z)r i - [P T (x) x 5(z-c)f n (5) 

After adding the contribution of steps located on the negative x axis, the interference function is obtained by Fourier 
transform of the auto-correlation function of the step positions : 

S(q x ,q z ) = FT{g(x,z)} 

= FT |«5(x) x 8{z) + J2 (PtW x S(z - c)]®" + [Pr(-x) x S(z + c)]®") | 
S(q x ,q z ) = l + (P T (q x )e^ c + (P T (q x )e^ c ) 2 ... + c.c.) (6) 



where Pr(qx) = J^°° PT(x)e l9xX dx is the characteristic function of the terrace widths distribution and c.c. is 
the complex conjugate part resulting from the contribution of steps located on the negative x axis. The geometric 
summation converges as | Priqx) |< 1 and one gets : 
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Where Re stands for the real part. Therefore the scattered intensity is : 
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This expression has been obtained by Pukite et al. (see Eg. 18-19 of Ref [15() or for instance by Rao et al. in the 



special case of the geometric distribution (see Eq. 16 of Ref [241]). Now an additional prefactor 



2sin(q~ c/2) 2 



(the CTR 



term) corresponding to the contribution of deep layers has to be included. Moreover the usual prefactor 4y is replaced 
by 



'-n 



to account for the discontinuous structure of the lattice. Equation [5] highlights the symmetric roles 



played by steps and terraces as e lqzC and Pr(qx) are respectively the characteristic functions of the step height and 
the terrace widths distribution. As shown in Appendix A, this result can be generalized for an arbitrary distribution 
of step heights Ps(z). 
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where Ps(q z ) is the characteristic function of the step height distribution. The pre-factors n ^„ }„ and 



2sin(q z c/2) allu 2 sin (<? x a/2) 

can be assigned respectively to the ideal truncation of an infinite crystal by a flat surface, the so called CTR (see Eq. 
[T]), and to the truncation of an infinitely large terrace by a step. In analogy to the CTR terminology, the 2s m(g a/2) 
term can be called Terrace Truncation Rod (TTR). Therefore, there is no intensity extinction at q x a = (2n+ l)n, i.e. 
in anti-bragg condition perpendicular to the step edges. Going further, a high sensitivity to the step edges roughness 
is expected by the combination of anti-bragg conditions both in q x and q z directions. This point will be discussed 
in more details within the two-level vicinal surface model. We have to add however that the scattered intensity may 
only be measurable for large enough miscut angles as the integrated intensity (i.e. in case of long range order) in 
these experimental conditions is proportional to the number of atoms located at step edges. 

To highlight the result given by Eq. let's consider the intensity scattered by a vicinal surface (lln) of a fee crystal 
(Fig[3]). The intensity is calculated in the plane of the reciprocal space perpendicular to the step edge direction (q x 0q z ). 
The atoms not located on the side view section (at y ^ 0) are included into the calculation thanks to the structure 
factor of the unit cell. Monoatomic steps are assumed with different terrace widths distributions. The distributions 
are chosen discrete by hypothesis to perform quantitative calculations in the whole reciprocal space (Appendix B). 
Of course if the atomic resolution is not required continuous distributions are enough. The additional shift in the 
x direction for adjacent terraces (5^75) is mandatory to keep the fee lattice structure. It is taken into account in 

the calculation adding a phase shift, e i9x 27I , to the characteristic function of the terrace width distribution. As the 
FWHM of the terrace size distribution gets larger, the scattering rod gets broader and second and third order satellites 
disappear in agreement with the increase of disorder. Of course depending on the miscut angle, the scattering rods 
are more or less tilted as they follow the direction perpendicular to the macroscopic surface plane. This remark is 
absolutely true close to a Bragg peak, however nonlinear rods may exist in some cases at larger wavevector transfer. 
For instance it will be shown at the end (Fig [15|) that a non integer value of the mean terrace size yields a S shape of 
the scattering rods. 



The two-level model 



The two-level model has been widely used in the literature to describe the intensity scattered by a surface during the 
homo-epitaxial growth in the sub-monolayer regime. The two level model allows evaluating the intensity scattered by 



a surface with a small roughness. The first models developed (e.g. (3 model [12L l25l|) do not provide any information 
on the line shape of the diffuse scattering parallel to the surface except a so called correlation length. Improvements 
have allowed analyzing the line shape assuming a terrace widths distribution 

0, 

16, 22], In this section, we revisit 

the previous calculations. The two-level model consists in alternating upward and downward steps (FigdJ). As for the 
vicinal surface, we consider the same elementary object, i.e. a semi-infinite sum of unit cells parallel to the terrace 
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plane. That way we can built step by step a two-level model folding nodes alternatively with this elementary object 
and its exact opposite (opposite electronic density). Their form factors are opposite of each other and reads: 



F±(q x , q z ) = ±F ceU (q x , q z ) ^ e mq * a = ±F ceU (q x , q z )- 



-iq x a/2 



n=0 
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(10) 



Assuming that P+(x) and P- (x) are respectively the upper and lower terrace widths distributions, the auto-correlation 
function of the electronic density reads for an upward step F + (x, z) as starting object and for x > 0: 



F + {x. z) ® F + (x,z) ® 5{x) + 

+ F + (x,z) ® F_{x,z) <8> P+{x) + 
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The prefactor 1/2 is used because the upward step contribution to the auto-correlation function as starting object 
has to be shared exactly with the downward step contribution. Therefore a similar g~(x, z) term is added considering 
a downward step F-(x, z) as starting object. Adding also the contribution of the steps located at x < 0, then the 
scattered intensity reads : 



I(q x ,q z ) = 2Rel± 
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A similar result has been obtained in the case of the stacking of layers [l!J, [26[ , and the exact result has already been 
found for instance by Pukite et al. (see Eq. 35 of Ref [IH), as well as by Croset and de Beauvais (see Eq. 20 of 
Ref [3). In these elder results, an additional [1 — cos(q z c)] dependence of the scattered intensity was obtained which 
originated from the hard-wall model hypothesis used for Helium scattering. For X-ray scattering it is exactly canceled 

yz, i.e. the structure factor of the semi- infinite column perpendicular to the surface 



by the CTR prefactor, 2sinfee/2) 2 
plane (Fig [5]). LEED patterns should be described as an intermediate case because a penetration length of a few 
atomic layers must be considered in the calculation. One can also see from Eq. [12] that the X-ray intensity is not zero 
at q z = 0. This is always the case for Helium or electron diffraction. Qualitatively one can understand the difference 
that way: for X-rays at q z — , i.e. in constructive interference conditions, one has to sum the electronic density of 
the atomic planes at different z coordinates deep into the bulk. As they have different coverage (for instance in case of 
a multi-level surface) each atomic plane contributes differently and the final sum does not give rise to a homogeneous 
atomic plane a priori. Therefore some intensity is expected in q z =0. For LEED or Helium scattering as only the last 
atomic layer contributes, the sum over z of the electronic density of each atomic plane gives a perfectly homogeneous 
atomic plane and therefore no off-specular intensity in q x . 
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The above calculation for X-rays is restricted to the off-specular intensity (q x ^ 0). The specular intensity reads : 



I(^Qz)=\F cell (0,q z ) | 2 



a iq z c/2 



2i sin (q z c/2) 



+ i 



(13) 



where 9 is the coverage. It only depends on the upper and lower terraces mean widths, which can be expressed in terms 
of the derivative at q = of the characteristic functions of the terraces widths distributions : 8 = P^(0)/(P^_(0) + 
PL(0)). Due to finite instrumental resolution and finite size effects, the separation of specular and off-specular 
contributions may be difficult to achieve. Measurements closed to anti-bragg conditions provide a powerful tool to 
solve this problem as the intensity is alternatively maximum off-specular and on the specular rod depending on the 
coverage. We refer to the work of Croset and de Beauvais [16| for a detailed study of this case. 

The N-level model 



The previous analysis is extended further with the N-level model. The N-level model or multi-level model is an 
extension of the two-level model for an arbitrary number of levels. Despite the simplicity of the present approach, the 
calculation of the scattered intensity for any N value remains difficult to obtain explicitly. For the sake of simplicity, 
only the final result for N=3 is given (Fig [5]): 



I(qx,q z ) 



Fcell (Qx , Qz 



2 sin(q x a/2)2 sin(q z c/2) 
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where P\(q x ), P%{q x ) and P 3 {q x ) are the characteristic functions of the terrace widths distributions of the first, 
second and third height levels. They are separated by monoatomic steps of height c and the probabilities for an 
upward or a downward step at the intermediate level are identical and equal to 1/2. This result has been obtained 
by Pukite et al. (see Eq. 40-44-45 of ref [HI) without the usual prefactor for X-rays siTl ^ q 1 c / 2 )i for the deep layers 



contribution and with the prefactor 



instead of ■ 



which now fully accounts for the atomic structure of the 



(<3xa) 2 sin((j x 

crystal. The analysis of Eq. [T4l clearly shows two distinct contributions respectively with pre-factors (1 — cos(q z c)) 
and (1 — cos(2(7 z c)), which can be assigned to the interferences between the waves scattered by terraces separated 
by respectively one and two height levels. As for the two-level model, an additional calculation yields the specular 
intensity. 



I(0,q z )=\F cell (0,q z ) | 2 



with 



a iq z c/2 



2i sin (q z c/2) 
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(15) 



(16) 



02 and 9 3 are the coverage of the second and third levels written in terms of the characteristic functions of the 
terrace width distributions for q x = 0. For larger N values the increase of complexity of the calculation makes the 
Phase Matrix Method [23[ more appropriate. 



The randomly stepped surface : the oo-level model 



The above models (2 and N-level models) are limited to the analysis of a small roughness i.e. bounded to a few 
atomic layers. They are not suitable to describe the intensity scattered by a rough surface with a divergent height 



7 



correlation function. An analytic method is proposed here (Fig [7]) to calculate the scattered intensity assuming a 
random distribution of monoatomic upward and downward steps (height c) separated by terraces having a width 
distribution Pt{x) [H, 27 1. Due to the divergence of the height correlation function, the whole scattered intensity is 
diffuse and no specular component has to be considered. 

Similarly to the 2-level model, downward and upward steps are two complementary objects, which in this case do 
not alternate regularly but randomly with a probability 1/2. The auto-correlation of the electronic density calculated 
for an upward step as starting object and for x > reads : 



^F+{x,z)®F + (x,z)® [S(z) x S(x) + ( <5(z) + 2 J(z C) ) x /VU) 



+ 



-S(z + c) + 5(z) - 5(z - c) + 5(z - 2c) 



x P T (x) ®P T (x) + ... 



^llib(l)(- S ( z + ( n 2fc ) c ) + S(z+(n-l- 2k)cj) I x P T (xf n + 



\k=0 



F(x,z)®F(x,z)® 5(z) x 5(x) + 



™ ELo 0) (-<?(* + (n 2k)c) + S(z + (n - 1 - 2fc)c)) 



n=0 



(17) 



where z) is the shape factor of an upward step, i.e. a semi-infinite row of atoms parallel to the terrace plane 

with positive electronic density. As for the two level model the downward steps is the opposite object, i.e. a semi- 
infinite row of atoms parallel to the terrace plane with negative electronic density F_(x, z) — —F + (x,z). Therefore 
downward steps contribute in the above expression of the autocorrelation function of the electronic density as negative 
terms. The expression of the auto-correlation function is not very intuitive but the advantage of this approach is to get 
ride of the terrace width distribution. (^J is a binomial coefficient which accounts for the all the possible combinations 
to reach a given height assuming random upward or downward steps. 

The expression of the intensity is obtained by Fourier Transform of the auto-correlation function of the electronic 
density after adding the contribution of downward steps as starting object (replacing S(z ± nc) by 5(z =F nc)) in the 
above expression and then including the terms for x < 0, i. e. the complex conjugate: 



I(q x ,q z ) = | F{q Xl q z ) | 2 x 
x U 



1 sin ^yY"^ Pr{qx) [ P t(<1x) cos (q z c) 
= \F(q x ,q z ) | 2 (l 



n=0 

ism(*f)e i3 ¥P T (q x ) 



+ c.c. 



1 - Pr{qx) cos (q z c) 



(18) 



In a more compact way: 



i(qx,q z ) 



I F ce ii(q x ,q z ) | 
2 sin (q x a/2)2 sin (q z c/2) 



Re < 2 



(1 - COS (q z c))(l - P T (q x )) 
1 - cos (q z c)P T (q x ) 



(19) 



This result is in agreement with the generalized equation of the scattered intensity by a vicinal surface given in Eq. 
[9l The cos (q z c) term in Eq. [19] is the characteristic function of the step heights distribution, which is simply in this 
case the sum of two delta functions in — c and +c weighted by 1/2. Therefore Eq. [HI which is given for a vicinal 
surface, can, in fact, be applied for any step height distributions including negative height, i.e. downward steps. This 
equation but without the s ^ q ] c / 2 )^ prefactor has been first obtained by Lu and Lagally (see Eq. 13 of Ref [28j) and 
then by Presicci and Lu in the special case of a geometric terrace width distribution (see Eq. 5 of Ref (2§|). Note that 
the scattered intensity in anti-Bragg condition (q z c — w) is the same for the randomly stepped surface and for the 
two- level model at half coverage as already mentioned by Henzler (30| . This ambiguity puts in evidence the necessity 
to measure the scattered intensity at different reciprocal space positions to fully characterize the surface roughness. 
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For instance, in Fig the q x position of the maxima corresponding to preferential terrace-terrace distances is a 
function of q z . This result can be compared successfully to Fig. 5 of Pimbley and Lu [31| calculated for a truncated 
geometric terrace size distribution. This can be understood through the interferences between the waves scattered by 
the terrace at different heights. This behavior does not occur for the two-level model for which the scattered intensity 
(off-specular) does not depend on q z . 

To characterize a rough surface, small angle scattering techniques are a common probe as they are preferentially 
sensitive to large length scales of the surface topography. In the small momentum transfer limit the intensity decreases 
as: 



I F ceU (q Xl q z ) 
{q x Df + 



(g*c) 4 



(20) 



where D is the mean distance between nearest neighbor terraces. This result can be compared with a more general 
approach based on the use of the height correlation function >, and in the special 

case of a ID system. 



F C eu(qx,qz) 

2 sin (q z c/2) 



-Hoc- 



-(l-co S (q z c))<(h(x)-h(0)) 2 > dx 



(21) 



To obtain an explicit expression the height correlation function has to be known. In the case of the paracrystal 
model, the asymptotic behavior of the height correlation function is a linear function of the distance if the first moments 



of the terrace width distribution exist (the distribution decays fast enough): < (h(x) — h(0)) 2 >- 



where c and 



D are respectively the mean step height and the mean distance between adjacent steps. In this case the calculation 
of the intensity from Eq. [HU reads exactly as Eq. [201 Of course the main drawback of the paracrystal model is its 
one-dimensional character contrary to the calculation based on the height correlation function that holds naturally 
to two dimensions and is more appropriate for instance to describe roughening 32, 3(|. However for some cases as 
proposed in paragraph 3, the paracrystal approach can be used to describe 2D systems. The present calculation 
has also the advantage to evaluate the scattered intensity in the whole reciprocal space (not only in the small angle 
regime) and for any terrace width distributions. 



The faceted surface 



Surface instabilities towards faceting or step bunching have been widely studied for instance in the case of the 
electromigration- induced faceting of silicon vicinal surfaces [37I ] or under adsorption of foreign species [38] . A model 
is developed here to calculate the intensity scattered by a faceted surface exposing two types of facets with two sets 
of size distribution. As an example, one of the facet type is a dense crystallographic plane and the other one is a 
vicinal face. To make simple calculations, the angle between the terrace plane of the vicinal facet and the low index 
crystallographic plane of the other facet is either 0° (Fig [9]) or 90° (Fig [101) an d respectively referred to small- angle and 
large-angle faceting. Similar calculations can be performed for any other angle. The principle of the calculation is to 
consider an elementary object which is a vicinal facet semi-infinitely extended in the x-direction (resp. z-direction) 
for the small-angle (resp. large-angle) faceting. For a vicinal (10N) facet with S steps, in the small-angle faceting 
case, the form factor reads: 



Ffacet(qx,qz) = F ce u(q x ,q z ) 
= F ce ii(q x ,q z ) 

And the interference function is (See Eq. [7]): 



t i7iq x a 



-iq^a/2 



S-l 



-2i sm(q x a/2) 



■ ( (q x N a+q z c)S \ 

Aqa:Na + qz c)(S-l) SUl(-^ '—) 

sin ( (gxAfa+ g . e ) y 



(22) 



S(q x ,q z ) = Re 



l + e i{l<cNa+q zC )Sp D ^ 
1 _ e i(q x Na+q z c)Sp D ( qx J 



(23) 
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Therefore the scattered intensity reads : 



F C eii(qx,qz) 



(2sin((faa/2)2sin( ( ^ a 2 + ^ c) )) 



Re ( 2 



1 _ e i(q*Na+q z c)Sp D ( qx ) 



(24) 



The previous result can be easily checked assuming that S = 1 and N = which gives exactly the intensity scattered 
by the vicinal surface (Eq. [H]). A more general expression of Eq. [M]can be obtained putting instead of e 4 ( 9xJVa+ ^ c ) s ', 
the characteristic function of the distribution of the number of steps, Pg, per vicinal facet. It results in the following 
expression of the intensity : 



F ce ii{q x ,q z ) | 



(2sin( ga fl/2)2sin( (g '* a 2 +9 - c) )) 



Re (2 



(1 - P S (q x Na + (fc C ))(l - Pp(q x )) 
l-P S (q x Na + q z c)P D (q x ) 



(25) 



For the large-angle faceting a very similar expression is obtained 



i(qx,qz) 



Fceii(q x ,q z ) I 2sm(q x Na/2) 



(2 sin (q x a/2)2 sin (q z c/2)2 sin( 



)) 



Re i 2 



(1 - P s {g x Na + g z c))(l - Pp(g z )) 
1- Ps(qxNa + q z c)P D (q z ) 



(26) 



In Fig [TTJ a simulation of the scattered intensity close to a Bragg peak shows the main features of the diffuse 
scattering by a faceted surface. The regularity of the faceting gives rise to satellite peaks nearby the Bragg peak 
whereas at larger wavevector transfer, scattering rods perpendicular to facets are clearly visible. 



TOWARDS TWO-DIMENSIONAL ANALYSIS OF SCATTERING FROM SURFACES 



The simplicity of the previous calculations arises because of the one-dimensional character of the surface morphology. 
However, this is also its limitation since in general surface topography should be described in two dimensions. The 
appropriate generalization of this method to two dimensions is the Markov Random Field theory as pointed out by 
Lent and Cohen [l4j . Without answering this question in the general case, a method is presented below to calculate 
the intensity scattered by 2D surfaces in few important cases : the kinked vicinal surface, the two-level vicinal surface 
and the step meandering on a vicinal surface. 



The kinked vicinal surface 



Analogously to the model of the vicinal surface the first step of the method is to find an elementary object that 
repeats regularly on the kinked vicinal surface. This object is a semi-infinite row parallel to the step (along the y-axis), 
starting at the kink (Fig [T2"]) . For instance, in the case of a monoatomic kink, the form factor of the elementary object 
is : 



o 

Pkinkcd vicinal Pee// (jlx i qy > ?2 ) ^ ^ e ^ y ^ 



n— — oo 



F ce u(q x ,q v ,q z )e l "y b / 2 

2ism(q y b/2) { ' 



In order to build the whole kinked vicinal surface, an infinite assembly of these objects are arranged along two 
main directions : the step edges (by definition) and perpendicularly to the steps to take into account any spatial 
correlations between kinks of neighboring steps. Both directions are decoupled by hypothesis and the kinks positions 
are given by the ideal 2D-paracrystal model (l9l. l39j|. assuming distance distributions between nearest neighbor kinks. 
However kinks (and steps) interact with each other via long range interactions (elastic, electrostatic) that yield a 
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well defined separation between neighboring kinks. To solve this problem it would be necessary to take into account 
correlations with far neighbors which would make the calculation much more difficult. Our purpose is to give a 
simple and analytical calculation of the intensity scattered by a kinked vicinal surface assuming generic distributions 
of distances between nearest neighbors kinks. Of course a relevant choice of the distribution laws may allow partially 
taking into account the kink interaction. 

In the direction parallel to the step edges the interference function is : 



S\\(qx,9y) = Re 



l + e^ a P Sl {q v ) 



(28) 



where Ps 1 {q y ) is the characteristic function of the distance separating two nearest neighbor kinks along the step 
edge and a is the kink width (monoatomic in this case) . This calculation is identical to the one for the vicinal surface 
model (see Eq. [7]). 

Similarly in the direction perpendicular to the steps, the interference function is given by: 



Sx(qx, q y , q z ) = Re 



\-e^P T {q x )P S2 {q y ) 



(29) 



where Pr(q x ) and Ps 2 (q y ) are respectively the characteristic functions of the terrace widths distribution and kink 
positions along the upper (lower) step edges. Ps 2 (q y ) denotes the variations of the kink positions along adjacent step 
edges. 

The expression of the intensity reads : 



I{q x ,q y ,q z ) = | linked vicinal xS\\(q x ,q y ) X Sj_(q x ,q y ,q z ) 

~ i 2 

F C eu{q x , q y ,qz) 



2sin(q y b/2) 



x Re 



l-e^PsM. 
l + e^PT{q x )Ps 2 {q y ) \ 
1 _ e^PMPsM J 

I Fceii (q xiq y -,qz) I 
2 sm(q x a/2)2 sm(q y b/2)2 sin(g z c/2) 



x fe 2 (l-e^)(l-P gl ( gg) ) 1 
i l-e^ a PsMy) J 
' (1 ~ e^ c )(l ~ P T (q x )P S2 (q y )) ' 
l-e^P T (q x )Ps 2 (q v ) 



x Re < 2 



(30) 



The first term in Eq. [30] consists of the multiplication of four terms : the form factor of the unit cell; the CTR, 
1/2 sin (q z c/2); the Terrace Truncation Rod, 1/2 sin (q y b/2), which has already been highlighted in the case of the 
vicinal surface model and a new term, 1/2 sin (q x a/2), which can be assigned to the truncation of an infinite step 
edge by a kink. Following previous terminology, it can be called Step Truncation Rod (STR). One should also notice 
that the scattered intensity is calculated in the whole reciprocal space (no specular scattering) as the whole sample 
volume (not only the surface) is taken into account in the calculation. As expected the intensity calculation shows 
two sets of scattering rods that are not perpendicular to each other and which correspond to the steps and kinks 
networks fFig [T3|) . Due to the miscut angle of the surface, the scattering rods shift along q z , in the q x direction (steps) 
as well as in q y direction (kinks) in order to be perpendicular to the average surface plane. The q y shift is hardly 
visible as the orientation of the normal to the macroscopic surface plane depends slowly on the kinks density. For 
completeness about Eq. 1301 a kink distribution can be considered (e.g. for non monoatomic kinks) replacing e %q * a 
by the characteristic function of the kink width distribution. The same can be performed for the step height (e ,<?zC ). 
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In addition to the topographic characterization of a kinked vicinal surface, this model can be applied to calculate the 
scattered intensity during epitaxial growth at kink edges. The form factor of the elementary object has to be modified 
according to : 



o 

iqr 



-^kinked vicinal — F ce i[{c[ x ; Qy ? Qz) ^ ^ 6 ^ y ^ -|- F a e 

n— — oo 

F ceU {q x ,q y ,q z )e^ 2 



2i sin(q y b/2) 



F a e« r (31) 



where F a is the structure factor of the additional atom located at a vector f from the kink edges. Therefore 
intensity variations might be observed at anti-Bragg conditions, i.e. in the center of the reciprocal lattice nodes 
{q x a — q y b — q z c — tt). However this may be hardly measurable due to the very small number of atoms involved in 
the scattering process resulting from destructive interferences between the waves scattered by all the atoms. Indeed 
without foreign atoms the intensity results from only half of the atoms localized at the kinks positions. This is similar 
to X-ray scattering by a flat surface in anti-Bragg conditions where only half of the atoms of the last surface plane 



give a contribution (cubic case, |34[). In homoepitaxy, no intensity change is expected (the kinks cannot be rough). 
However in heteroepitaxy, the intensity can vary significantly depending on the structure factor of the foreign atom 
at its position with respect to the kink site. 



The two-level vicinal surface 



The aim of this part is to calculate the intensity scattered by a vicinal surface which has rough step edges [fj, |4fJ • 
For that purpose we consider that the step edges can be described by a two-level model (Fig PH)) . Making use of 
previous results for the two-level model and for the vicinal surface the scattered intensity reads : 



i(q x ,q y + o,q z ) 



x Re 



x Re 



Fceii(qx,q v ,qz) 

2sm(q y b/2) 

' (l-P + (g,))(l-P_( g ,)) 
1 - P+(q y )P-(q y ) 

l + e^ c P T (q x )P S2 (q y ) \ 
1 - eJ^P T (q x )P S2 (q y ) J 



(32) 



Where the first two terms at the right hand side of the equation has already been given in the case of the two-level 
model (see Eq. |T2"j) . Indeed a two-level model for step edges is equivalent to a one dimensional two-level model for a 
surface. The third term deals with the interference function for adjacent steps (see for instance Eq. [7] of the vicinal 
surface) with an additional term Ps 2 (ly) which accounts for correlations between the positions of kinks at adjacent 
steps. 

In this case the calculation of the scattered intensity is restricted to q y ^ 0. Similarly to the two-level model one 
has to make an additional calculation for this case : 



I(Qx,0,qz) =| F ceU (q Xl 0, q z ) 



-iq m a/2 



-2i sm(q x a/2) 



+ 9e~ tq * a 



Re 



1 + e^ c P T (q x ) 
1 _ e «fecp T ( gx ) 



(33) 



Where 9 is the analogue of the coverage applied to a step edge. As already mentioned in the part dedicated to the 
one-dimensional analysis of vicinal surfaces, the roughness on step edges modifies the Terrace Truncation Rod in a 
way similar to the roughness on a flat surface for the CTR (Eq. [33l) . This remark calls for additional comments. For 
instance in the case of homoepitaxial growth on a vicinal surface, intensity oscillations may be observed in anti-Bragg 
conditions if the step edges are alternatively rough and smooth (Tig |15|) . To maximize the sensitivity to the roughness 
one has to be in anti-Bragg condition on the scattering rod of the step network (q x a = n, q y b — 0, q z c = n). Fig 
[T51 reveals also that the scattering rods have a S shape. This can be assigned to the non integer value of the mean 
terrace size. Indeed, close to a Bragg peak, the scattering rods are perpendicular to the average surface plane. As 
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the mean terrace size is not an integer, the scattering rods do not point towards a Bragg peak. Therefore there is a 



slight orientation change close to anti-Bragg conditions to connect two head-to-tail scattering rods [31(. In addition 
shoulders in q y direction originating from a preferential distance between kinks are expected in the diffuse scattering 
(Eq. 152")) as pointed out by Wollschlager et al. and revealed experimentally by Tegenkamp et al. [|| with SPA- 
LEED. This effect should be enhanced at half coverage of the step edges (maximum roughness) exactly in the same 
way as the two-level model for a surface (see Ref [Hj]). In the case of heteroepitaxial growth at step edges on vicinal 
surfaces (e.g. Co/Pt(997) [4l[) complex interference effects are expected arising from the different scattering power 
between the substrate atoms and the adatoms. This model can be extended to the N-level vicinal surface. One has to 
add that the correlation between the morphologies of neighboring step edges may be adjusted thanks to the Ps 2 {<ly) 
function. It is related to the kink-kink relative displacement along_ two adjacent step edges. If the kinks positions 
are fully correlated, e.g. at identical y position from step to step, Ps 2 { c Ly) — 1- If the morphologies are independent 
Ps 2 ( ( ly) — Hly)- 111 this latter case, for q y ^ the scattered intensity is the incoherent sum of the contributions 
of independent step edges. For a detailed study of this case we refer to the work of Wollschlager et al. [40]. Any 
intermediate case can be considered. 



Step meandering on vicinal surfaces : the oo-level vicinal surface model 



The last example deals with step meandering on vicinal surfaces [32|. Instabilities on vicinal surfaces induced 
by specific growth conditions or an external stress may lead to step meandering [42l |. Making use of previous re- 
sults on randomly stepped surfaces and vicinal surfaces, the generalization to step meandering on vicinal surfaces is 
straightforward (Fig [T6|) . The scattered intensity for a kink size a and a step height c is : 



I F ceU (q x ,q yi q z ) \ 
2 sm(q x a/2)2 sin(q y b/2) 



x fle 2 (l-co S (M)(l r ^fa)) >x 



x Re 



1 - cos (q x a)P Sl (q y ) } 

(1 - e^ c )(l - P T (q z )PsM) 
1 - e^P T {q x )PsMy) 



(34) 



Where the first two terms on the right hand side of the equation are directly extracted from Eq. [19] of the randomly 
stepped surface. Indeed a step edge with random kinks is equivalent to the one dimensional randomly stepped surface 
model. The last term is the interference function of adjacent steps as seen previously for instance in the case of the 
kinked vicinal surface. A simulation of the intensity scattered by the step meandering model is shown in Fig 1 171 The 
intensity decreases much faster in the q x direction compared to the two-level vicinal surface model (q x instead of 
q x 2 ) due to very rough step edges. 

For this model one may wonder if steps may collide which is a non physical situation due to step-step repulsion. 
Even if the step distance distribution, i.e. Pt(x), takes only positive values, steps may collide due to the Ps 2 {y) 
probability distribution for the kink-kink distance of adjacent steps. The effect of Ps 2 (y) is to decorrelate the kink 
positions of two neighbor steps. If one choose a very broad Ps 2 (y) probability distribution, a step may collide with 
an adjacent step. However this can be avoided by choosing a Ps 2 (y) probability distribution with a maximum | y \ 
displacement much smaller than the product of the kink-kink distance on a step edge times the number of kinks to 
fill the terrace width. This condition is not very restrictive but one can see that the limit case Ps 2 (y) = 1 does not 
fulfill this requirement. It leads to the diffraction by decorrelated steps as pointed out by Wollschlager et al. [ijj], but 
of course uncorrelated steps with random kinks on the step edges may collide. This remark on P$ 2 (y) applies also 
for the kinked vicinal surface model but not for the two-level vicinal surface model as the step-step fluctuations are 
confined in ±1 atomic distance and the non crossing condition for steps is easily fulfilled. 

As for the kinked vicinal surface, the calculation provides the whole scattered intensity as their is no specular 
scattering. Of course the above expression can be generalized for any kink size and step height distributions. For 
instance it is straightforward to combine the randomly stepped model with a random distribution of kinks along the 
step edges just changing e lqzC by cos (q z c) in Eq. [MJ 
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LATTICE DISTORTION 



Until now no lattice distortion have been introduced and the crystallographic structure is the same at the surface 
and in the bulk. This assumption may be the most limiting one to quantitatively characterize surfaces with X-rays 
as steps and kinks are intrinsic surface defects which stress the sample. To solve this problem, we suggest 



an approximate calculation of the scattered intensity which takes into account the strain field induced by steps and 
kinks assuming that the surface morphology is not perfectly periodic. To illustrate this point, let's consider the case 
of a vicinal surface. For a perfectly periodic network of steps the strain field along only one atomic plane parallel 
to the terrace and starting at a step edge allows describing the strain field inside the whole sample volume [Zij ]. 
This aforementioned atomic plane is nothing else than the semi-infinite elementary object used to build the vicinal 
surface model (see Fig [2|. Consequently the strain field calculation (e.g. in the framework of the linear elasticity 
theory) obtained for a perfect network of steps, may be applied to a slightly disordered one by simply using the form 
factor of the strained elementary object. Of course for highly disordered vicinal surfaces this calculation will not be 
appropriate. One has to be aware that a variation of the terrace size modifies not only the amplitude of the strain 
field but also the positions of the accidents of the displacement field which are located below the nearest neighbor 
steps. Therefore much caution is necessary to use this approximation. A similar result can be obtained for the kinked 
vicinal surface. The strain field induced by kinks on an ideal kinked vicinal surface (bi-periodic kinks network) can be 
obtained through the strain field calculation along only one atomic row parallel to a step edge and starting at a kink. 
Therefore an approximate calculation of the scattered intensity by a strained kinked vicinal surface can be obtained 
transferring this latter result to the elementary object. However contrary to the vicinal surface case, the strain field 
induced by an ideal kinked vicinal surface has not been calculated yet. 

CONCLUSION 

We have developed a general simple method to calculate the intensity scattered by different kinds of surface 
morphologies that are not perfectly ordered. The model based on the paracrystal framework allows to recover 
previously published results on the scattering from vicinal surfaces, two-level models and randomly stepped surfaces. 
The most important result is the ability to characterize two dimensional surface structures from nearest neighbor 
distance distribution between steps or between kinks. The main hypothesis of the model consists in neglecting the 
long-range correlations between neighboring steps (or kinks) and, for the two dimensional models, in assuming that 
the step edge shape and the terrace width are independent. In the future, additional correlations should be considered, 
e.g. in the field of thin film growth between the size of neighboring terraces as highlighted first by M. C. Bartelt and 



J. W. Evans [47|. A step in this direction has been proposed in the case of the two- level model (see Ref [48[) which 



should be extended to more complicated surface morphologies. 
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APPENDIX A: INTENSITY SCATTERED BY A VICINAL SURFACE WITH A STEP HEIGHT 

DISTRIBUTION 

In this section, the intensity scattered by a vicinal surface with both step height Pjj and terrace width Pjj distri- 
butions is calculated (Fig [T8|) . 

The considered scattering object is a packing of H crystallographic planes, semi-infinite and parallel to the terrace. 
H is distributed according to the probability law P#. 
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The auto-correlation function of the vicinal surface is (only the contribution of steps at x > is given): 
g+(x,z) = JJJJ [F (x,z)® F (x,z)®d{z)®d{x) + 

+ F Q (x,z)®F 1 {x,z)®8(z °— — -)®P Da {x) + 

+ F (x, z) ® F 2 (x, z) ® S(z - H ° + 2 H2 ) <8> P Do (s) <8> (x) ® P Dl (*)] + ... 

+ F (x, z) ® F n (ar, z) ® <5(z - H °+ Hn ) ® p Dq ( x ) ^ [p H ( x ) ^ P D (x)]»" _1 + ... 
+ . . . ] dP Do . . . dP Dn _ , dP Ho . . . dP Hn . . . 



(35) 



where F n (x, z) is the shape of the n th object, i.e. F n (x, z) = 1 at the positions of the unit cells otherwise it is zero. 
The intensity is obtained by Fourier transform of the auto-correlation function : 



I(q x ,q z ) =FT{g(x,z)} 
Fo I 2 + 



H 

+ 
+ 
...+ 

...+ C.C. 



F Q {q x ,q z )F*(q x , qz y^ Jl£i ^^^ D ° a + 



Fo(q x ,q z )F;(q x ,q z )e^^^^e^ D " a e^ D ^ H ^ + ... 



H ,H 1 ,H 2 ,D ,D 1 



II,. // ,D /)„ : 



dP Do ...dP Dn _ 1 dP H() ...dP Hn ... 

As the structure factor of the object is: 

F n (q x ,q z ) = F ceU {q x ,q z )- 
= F ceU (q x ,q z )- 



(36) 



sin (q z H n c/2) 



-2isa\{q x a/2) sin(g z c/2) 

e -iq x a/2 e l(q z H 7l c/2) _ e -i{q z H n c/2) 

-2ism(q x a/2) 2isin(g z c/2) 



(37) 



Therefore the scattered intensity reads: 



Fceii(q x , q z 



(2 sin ((^0/2)2 sin (g z c/2)) 

(PH{q z c)-l)(l-P H (q z c))P D {q x a) 



2Re y. - P H {q z c) + 

Fceii{q x ,q z ) I 



1 - PH(q z c)P D {q x a) 



(2 sin (q x a/2)2siri {q z c/2)) 



Re {2 



(1-P g fec))(l-Pr>(g*a)) 

l - PH{q z c)P D {q x a) 



(38) 



which is the expected result. 



APPENDIX B: DISCRETE DISTRIBUTIONS AND CHARACTERISTIC FUNCTIONS 



Some examples of discrete distributions are given as well as their characteristic functions 
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The Binomial distribution 



The Poisson distribution 



The Geometric distribution 



PpA") = c n NP n (i-p) N ~' 
p p , N (q) = {i- P + ve iq ) N 



Pu{n) = - 
n 



Pu{q) 



\ 



P p (n) = p(l-p) n 
Pp{q) = l-(l-p)e"i 

For completeness, more complex distributions can be considered such as the hypergeometric distribution or the 
negative binomial distribution. 
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FIG. 1: The vicinal surface described as a stacking of terraces. Pt(x) is the terrace width distribution and c is the step height. 
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FIG. 2: The vicinal surface model, a) The elementary object is a semi-infinite plane starting at the step edge, b) and c), 
the vicinal surface is built step by step assuming a terrace width distribution Pt(x). D is the average terrace width and a is 
the FWHM. The probability laws of having a distance separating the first, second and third neighbor steps are schematically 
drawn on the x axis. Due to an increase of disorder the probability laws broaden with the distance. 
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FIG. 3: (a) Side view of a (115) vicinal surface of a fee crystal and scattered intensity (logarithmic scale) in the (q x 0q z ) plane 
for different terrace width distributions. From (b) to (d) the terrace width distribution gets larger and larger, making the 
scattering rods less and less narrow, a is the lattice parameter, a/2 is the step height and a/(2\/2) is the shift in x of the lattice 
parameters of two adjacent terraces. 
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FIG. 4: The two-level model is built step by step assuming alternating upward and downward steps. Terrace width distributions 
are respectively P+(x) and P-(x) for the high and the low levels. D± and a± are respectively the average terrace width and 
the FWHM for both distributions. The probabilities of finding a distance x separating a step up from a step down (or up) are 
schematically drawn on the x-axis. 
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b) 



FIG. 5: a) Hard- Wall model (e.g. used for Helium Atom Scattering) for which only the last atomic layer is taken into account 
in the calculation of the scattered intensity. It results in an additional (1 — cos(g z c)) term in the calculation of the scattered 
intensity, b) Model for X-ray scattering for which the whole sample volume has to be considered due to the very weak cross 
section. The difference in the off-specular scattering calculation consists in considering the atomic layer below the upper 
terraces. 
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FIG. 6: The 3-level model is built assuming that from level 1 and level 3, the neighboring steps must be respectively upward 
and downward. From level 2, both upward and downward steps are equally probable. The terrace width distributions are 
Pi (a;), P2{x) and Ps(x) for each level. 
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FIG. 7: The rough surface can be seen as an infinite level model. The terrace width distribution is Pt{x) D is the mean distance 
and a is the FWHM, and c is the step height. Upward and downward steps are equally probable. As for the vicinal surface 
and the two-level model, the distributions of distances separating a n th neighbor step from an arbitrary step are schematically 
drawn on the x-axis. The FWHM of the distributions increases with the distance. 




FIG. 8: Reciprocal space map (q x ,q z ) of a rough surface. D is the average distance between terraces. a/D = 0.3 (binomial 
distribution, see Appendix B) 
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FIG. 9: Model of a faceted surface (small-angle faceting) showing two types of facets, a vicinal one of (ION) type and a low 
index one of (001) type. Two sets of size distributions are included: the number of steps Ps of the vicinal facet and the number 
of unit cells Pd of the low index facet. 
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FIG. 10: Surface model for a large-angle faceting. Same as before but the low index facet is a (100) type. 
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FIG. 11: Scattered intensity map (q x ,q z ) in logarithmic scale (a decade of intensity separates each color change) of a faceted 
surface (large-angle faceting type). q x is parallel to the macroscopic surface plane and q z is perpendicular. The vicinal facet 
is of type (102). A schematic drawing of the surface morphology and the size distributions of P(S) and P(D) are given. A is 
the average periodicity of the faceting and a is the angle between the macroscopic surface plane and the most extended facet 
(vicinal facet). 
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FIG. 12: The kinked vicinal surface model is built assuming that the elementary object is a semi-infinite row starting at the 
kink edge. Then the step edge is built from the kink- kink distance distribution, Ps 1 (y), and the terrace width is controlled 
by Pt(x). Ps 2 (y) is an additional distribution to control the correlations of kinks positions of neighboring steps, a and c are 
respectively the kink size and the step height. 
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FIG. 13: a) Step distance distribution Pt(x), kink-kink distance distribution along the step edge Ps 1 (y), and kink-kink distance 
distribution along nearest neighbor steps Ps 2 (y) ■ b) Simulation of the corresponding surface morphology, c) Intensity map 
scattered by this kinked vicinal surface at q z — 0. 




FIG. 14: The two-level vicinal surface model allows describing the step-edge of a vicinal surface as a two-level model. P+(y) 
and P-(y) are the width distributions of the higher and lower levels. As for the kinked vicinal surface model, the terrace width 
is controlled by Pt{x) and the kink-kink distance along the y-direction for neighboring steps by Ps 2 (y). 
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FIG. 15: a) Intensity map (logarithmic scale) of the two-level vicinal surface model in the (q x 0q z ) plane, b) Histogram of the 
terrace width used for the calculation, c) Intensity oscillations calculated in anti-Bragg condition (white cross on a). 9 is the 
coverage rate of the step edge. One can notice the S shape of the scattering rods arising from the non integer value of the mean 
terrace size. 




FIG. 16: The step meandering model on a vicinal surface is built assuming that the step edge can be described by the infinite 
level model (kink-kink distance distribution: Ps 1 (y)). The terrace width is controlled by Pt{x) and the kink-kink distance 
along the neighboring step edge in the y-direction by Ps 2 (y). 
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FIG. 17: a) Intensity map (logarithmic scale) ol the step meandering vicinal surface model in the {q x Qq z ) plane, b) Histogram 
of the terrace width distribution, c) Histogram of the kink- kink distance distribution : Ps 1 (y). Ps 2 (y) is identical to Ps 1 (y) 
but centered around y = 0. Compared to the two-level vicinal surface, the intensity decreases much faster in q x because the 
step edges are much more rough. 
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FIG. 18: Vicinal surface model assuming a terrace width distribution (Pt) as well as a step height distribution (Ph)- D and 
H are the mean terrace width and the mean step height and gt and oh are the FWHM of the probability laws. The FWHM 
of the probability to find two steps separated by a certain distance increases with this distance as schematically illustrated on 
the x-axis (same for the height, z-axis). 



